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ABSTRACT 

Infrared imaging properties of dusty winds around late-type stars are investigated in 
detail, employing a self-consistent model that couples the equations of motion and 
radiative transfer. Because of general scaling properties, the angular profiles of surface 
brightness are self-similar. In any given star, the profile shape is determined essentially 
by overall optical depth at each wavelength and it is self-similarly scaled by the size 
of the dust condensation zone. We find that the mid-IR is the best wavelength range 
to measure directly the angular size of this zone, and from IRAS data we identify the 
15 best candidates for such future observations. We also show that the visibility func- 
tion at short wavelengths (5s 2 |a.m) directly determines the scattering optical depth, 
and produce theoretical visibility curves for various characteristic wavelengths and 
the entire parameter range relevant to late-type stars. The infrared emission should 
display time variability because of cyclical changes in overall optical depth, refiect- 
ing luminosity- induced movement of the dust condensation point. Calculations of the 
wavelength dependence of photometric amplitudes and time variability of envelope 
sizes are in agreement with observations; envelopes are bigger and bluer at maximum 
light. 
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1 INTRODUCTION 

Ongoing progress in imaging techniques, expected to con- 
tinue at an accelerated pace (e.g. Gezari & Backman 1993; 
Skinner et al. 1994), enables spatial resolution of infrared 
sources. The wealth of details provided by high-resolution 
data requires a commensurate theoretical effort, which can 
be expected in general to involve a large number of free pa- 
rameters. The freedom afforded by many parameters usually 
produces good fits but casts doubt on the uniqueness of such 
models. 

Late-type stars are a notable exception. Recently we 
have shown that, given the grain composition, the IR signa- 
ture of dusty envelopes around such stars is essentially de- 
termined by a single parameter - the flux-averaged optical 
depth Tp (Ivezic & Elitzur 1995; hereafter IE95). This pa- 
rameter alone controls both the dynamics and IR emission. 
The flue resolution achieved in recent interferometric obser- 
vations (envelopes are resolved on angular scales as small 
as 0.03 arcsec at ~ 10 \Lm; see Danchi et al. 1994) mer- 
its detailed modelling of these stars, which we offer here. 
The surface brightness distribution can be calculated once 
the density and temperature distributions of the dust in the 
envelope are known (e.g. Crabtree & Martin 1979). In pre- 
vious studies, both distributions were usually modelled by 
power laws prescribed beforehand. Instead, here we deter- 



mine these distributions self-consistently from a complete, 
detailed model that solves the envelope dynamics and radia- 
tive transfer simultaneously. This self-consistent approach 
enables us to compare different types of observations and to 
investigate systematically which techniques are best suited 
for determining various quantities. 

General analysis of surface brightness distributions is 
described in Section 2. The effects of stellar variability on 
observed IR emission are discussed in Section 3. We analyse 
the observational implications of our results in Section 4 and 
summarize them in Section 5. 



2 SURFACE BRIGHTNESS DISTRIBUTION 

The surface brightness of a spherically symmetric source de- 
pends only on the impact parameter b, the displacement 
from the centre of symmetry. In the case of late-type stars, 
we have shown (IE95) that the solution of the radiative 
transfer problem is a self-similar function of the scaled vari- 
able b/ri, where ri is the envelope's inner boundary, corre- 
sponding to the dust condensation radius. For a given type 
of dust grains, this function is virtually fully determined by 
the overall optical depth r and does not depend separately 
on luminosity, mass-loss rate, etc. 

Fig. 1 presents the results of our model calculations at 
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Figure 1. Normalized surface brightness distributions at 2.2, 5 
and 10 i-im for models with dust condensation temperature T\ 
= 800 K, stellar temperature = 2500 K and optical depths as 
marked. Here h and 9 are, respectively, the linear and angular dis- 
placements from the centre, ri is the radius and 9i is the angular 
diameter of the dust condensation zone. 



the three representative wavelengths: 2.2, 5 and 10 \i.m. The 
surface brightness profile depends only on the dimensionless 
variable 6/ri, so that it scales self-similarly with the dust 
condensation radius r\. At any given wavelength, the shape 
of this profile is determined by r; two different systems with 
the same r will have identical profiles in terms of b/r\, so 
that the actual value of ri is irrelevant. For a system at 
distance D, the impact parameter b corresponds to an an- 
gular displacement 9 — b/D. If we introduce the angular 
diameter of the dust condensation zone Oi — 2r\/D, then 
the observed surface brightness is a self-similar function of 
fo/n = 29 /Ox. 

In these calculations we have used a detailed model that 
takes into account the wind acceleration by radiation pres- 
sure, determined self-consistently from a solution of the cou- 
pled equations of hydrodynamics and radiative transfer (for 
details see Netzer & Elitzur 1993; IE95). Thanks to scaling, 
the displayed profiles are applicable to any late-type star 
at the optical depths indicated for each wavelength, with a 
given star corresponding to different optical depths at dif- 
ferent wavelengths. The models were calculated with the 
opacities of amorphous carbon (Hanner 1988). Other types 
of grains (silicate, SiC, etc.) produce almost identical pro- 
files for the same optical depths. For a given star, only the 
association of optical depths with wavelengths may depend 
on grain chemistry when there are spectral features in the 
opacity. For example, for the same optical depth at 2.2 \im, 
a star with silicate grains will have a larger 10-p.m optical 
depth than one with carbon grains. Thus the two stars will 
have the same profile at 2.2 jim but different ones at 10 \im. 
The 10- [im profile will be the same for the star with sili- 
cate grains and another star with carbon grains whose r is 
matched at that wavelength. 

In addition to chemical composition, optical proper- 
ties depend also on the grain size a. Except for spectral 
features, the absorption and scattering efficiencies are self- 
similar functions of the scaled variable x — 2na/X and their 
functional form does not change when x < 1. Therefore, 
so long as only a; < 1 is covered by observations, the grain 
size is largely irrelevant. For the wavelengths addressed here, 
A ^ 2.2 \Lm, all grains smaller than ~ 0.2 \Lm are practically 
equivalent, and following IE95 we adopt a single grain size of 
0.05 [im. The possible effect of grains larger then 0.2 (j.m is 
discussed below (Section 2.1). Other input parameters used 
in these calculations are the stellar temperature T, — 2500 
K and the dust condensation temperature Ti — 800 K. The 
model results depend primarily only on the ratio of these 
two temperatures. Varying Ti and T* within their relevant 
ranges has a negligible effect on the results. 

The narrow central component evident in the profiles is 
the stellar radiation attenuated by the dust shell. Its width 
is r,/ri, where r* is the stellar radius. This ratio is directly 
related to T^/Ti. The results of our detailed model calcula- 
tions, which account for the contributions of both the central 



star and the dusty envelope to the radiation field at ri, can 
be adequately summarized with 

Here a is a dimensionless coefficient that depends primarily 
on the grains' optical properties and only weakly on the 
overall optical depth. For carbonaceous grains a is ~ 1.2 for 
T — 0, increasing slowly with r to a ~ 1.4 at r(2.2 \un) 50, 
the largest plausible optical depth. For silicate grains, the 
corresponding range is ~ 0.5-0.8. In addition, a has an even 
weaker dependence on Ti and T* , which can be neglected in 
all practical applications. 

The extended component visible in the figure is the con- 
tribution of the dust shell. At small r this intensity compo- 
nent is proportional to the optical depth along the line of 
sight, and the displayed profiles reflect the behaviour of this 
function. The features evident at fe/ri — ±1 simply reflect 
the envelope's inner edge. High-resolution imaging will de- 
lineate these features as a bright rim whose angular size (^i) 
is wavelength independent. The angular size of the dust for- 
mation zone is then directly measurable. In addition, the 
optical depth can be estimated by comparing the measured 
ratio of peak intensities of the central (stellar) component 
and the dust formation features with those of the displayed 
proflles. As the optical depth increases, these features grad- 
ually disappear and the shell's contribution approaches that 
of the star. The intensity profile becomes a featureless bell- 
shaped curve, whose width is roughly proportional to t9i. 
In this case it is impossible to determine r and 9i indepen- 
dently. 



2.1 Visibility curves 

Most spatially resolved IR observations of late-type stars 
are interferometric and produce the two-dimensional Fourier 
transform of the surface brightness, the visibility V{q), 
where q is the spatial frequency. For a circularly symmet- 
ric intensity profile, V{q) = F{q) /F{0), where 

oo 

F{q) = 2n J I{9)Jo{2nq9)9d9 (2) 



(cf. Rogers & Martin 1984; Jo is the Bessel function of zcroth 
order). In the present case 

V{q) = fsVs + .hVo, (3) 

where Vs and Vo are the individual visibility curves of the 
stellar and dust emission, respectively, and fs and /d are the 
fractional contributions of these two components to the total 
flux. The stellar component can be closely approximated by 
a uniform-brightness disc 

Vs = 'Jl^, (4) 
nq9t 

where Ji is the Bessel function of first order and 9t = 2r,/_D 
is the angular diameter of the stellar disc. The stellar visi- 
bility is a function of q9, . Because of scaling, the shell's con- 
tribution to the intensity is a self-similar function of ^/^i, 
thus the corresponding visibility Vy) is a self-similar function 
of q9x. 
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Figure 2. Visibility curves {q is the spatial frequency) at 2.2, 5 
and 10 \xm for models with the same temperatures as in Fig. 1 
and optical depths as marked. The dashed line is the visibility 
function of the naked star (r = 0). 



Fig. 2 presents the visibility curves of the surface bright- 
ness distributions displayed in Fig. 1, including more values 
of T for optically thin envelopes. Just as the surface bright- 
ness profiles scale self-similarly with ri, the visibility pro- 
files scale self-similarly with Oi. The dashed lines are the 
visibility curves of a naked star, Vs, included for compar- 
ison. They are quite flat, reaching zero visibility only at 
qOi = 1.2 01/6, = 1.2ri/r,. Note that the ratio n/r., typ- 
ically ~ 5-10, is directly related to the temperature ratio 
Ti/T. (equation |l|). 

The stellar and envelope contributions are easily dis- 
cernible in the figure. All visibility curves display an initial 
drop with qOi, reflecting the extended nature of the dust 
component, followed by a levelling off when the envelope 
is fully resolved and the stellar component takes over. The 
value Vc of the visibility at that levelling-off is roughly the 
fraction of the stellar contribution to the total flux. The sig- 
niflcance of this contribution depends on both optical depth 
and wavelength. When it is appreciable, K is flnite and the 
levelling-off occurs at qOi ~ 1. The visibility curve can then 
be used to determine both 9i and r - the magnitude of Vc 
directly gives r, and the location of the levelling-off deter- 
mines the angular scale 9i ~ 1/q. In sources where Vc — 0, 
the levelling-off occurs at progressively smaller qOi as r in- 
creases, a direct consequence of the broadening of the in- 
tensity profile (cf. Fig. 1). The visibility curves become self- 
similar with angular size-scale roughly proportional to t9i, 
thus T and di cannot be determined independently. These 
points have been recognized previously in the discussion of 
surface brightness. 

It is important to note that the dust component of 
the visibility profiles always varies with wavelength, even 
in optically thin envelopes whose images are dominated by 
wavelength-independent bright rims (corresponding to the 
dust formation zone). The reason is that, at any given wave- 
length, the visibility refiects the angular size of the region 
from which most of the fiux is emitted. This size corresponds 
to the radius where the dust temperature is comparable to 
the Wien temperature of that wavelength. Since the dust 
temperature decreases with distance, the size of the fiux 
emitting region increases with wavelength. Such variation 
has been detected in many observations. Most visibility ob- 
servations are interpreted in terms of Gaussian fits and are 
reported in terms of the resulting Gaussian FWHM 9g- A 
quantitative comparison of theory with observations must 
incorporate this widespread method of reporting observa- 
tional results. The observed Gaussian widths are found to 
increase with wavelength according to Gq oc A'', and the in- 
dex k is found to vary considerably among individual stars. 
In the wavelength range 2-10 \im, Dyck et al. (1984) ob- 
tained k ~ 0.2 for C stars (carbon dust) and k ~ 0.7-1.3 for 
O stars (silicate dust). Cobb & Fix (1987) found ~ 1 in 
the wavelength range 5-10 ^.m for a sample of OH/IR stars 
(silicate dust). 

Our models readily explain these detailed observations. 
Fitting a Gaussian to the model visibilities, we find that in 



the wavelength range 2.2-10 ^.m these fits can be summa- 
rized with 



A 



3.2 li-m 



r < 1 

T>1, 



(5) 



an approximation accurate to within 10-20 per cent. There- 
fore, in optically thin envelopes k = 0.6 independent of grain 
chemical composition, in agreement with the lower bound of 
the range found in O stars by Dyck et al. In optically thick 
envelopes, on the other hand, k does vary with the grain 
chemistry because of its dependence on t. We find that 
oa A°'^ in amorphous carbon optically thick envelopes, 
in agreement with the observations of C stars. The results 
for silicate grain optically thick envelopes are more involved 
because of the peak at 9.7 \im. A parametrization in the 
form 6g oc A'' results in fc = 1.3 when the fit applies to the 
wavelength range 2-10 [im and k — 4 for the range 7-10 
[im. These results are in agreement with the various O star 
observations quoted above. 

In addition, observations by Fix & Cobb (1988) at 8 
and 10 nm show that is roughly the same at both wave- 
lengths in IRC-l-10420 but increases by about a factor of 2.5 
in OH 26.5-f 0.6, even though both stars display the silicate 
feature. Our results explain both observations if the dust 
shell is optically thin in IRC-l-10420 (expected 6'g variation 
of only 10 per cent) and thick in OH 26.5-1-0.6 (expected 
variation by a factor of 2.4). Indeed, the former displays the 
silicate feature in emission, the latter in absorption, indicat- 
ing optically thin and thick envelopes, respectively. 

In spite of these successes, our model results indicate 
that single- Gaussian fits to visibility curves must be used 
with great care. Such fits are adequate within observational 
uncertainties (~ 0.05 in visibility) only at short wavelengths 
(Si 2.2 nm). The quality of Gaussian fitting deteriorates as 
the wavelength increases, and we find such fitting meaning- 
less for A J; 10 [im. The apparent success of Gaussian fits at 
such wavelengths in a number of actual observations usually 
refiects partial coverage of the visibility curve (e.g. Sutton, 
Betz & Storey 1979, V < 0.5; Cobb & Fix 1987, V > 0.8). 
Indeed, when the spatial frequency coverage provides the en- 
tire visibility curve, single-Gaussian fits become inadequate 
(e.g. Mariotti et al. 1983; Dyck et al. 1987). While such a 
fitting procedure fails, our models provide satisfactory fits 
to all observations and provide a more adequate method for 
data analysis. 

Finally, for given a grain size a, scattering is important 
only for wavelengths shorter than Asca ~ 2na. For the size 
used here of 0.05 |J.m, Asca is 0.3 \im and scattering is negligi- 
ble at all the displayed wavelengths. For any plausible grain 
size (Ss 1 vim), scattering can become important only at the 
shortest displayed wavelength of 2.2 ^.m. When that hap- 
pens, the shape of the visibility curves remains essentially 
the same but the corresponding optical depth is significantly 
smaller. For example, if the grain size is 0.2 \i.m, the five 
curves in the topmost panel of Fig. 2 are produced with op- 
tical depths of 0.1, 0.8, 1, 3 and 10, respectively, instead of 
the indicated values of 0.1, 1, 3, 10 and 30. 
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3 TIME VARIABILITY 

Late-type stars display variability with periods of order 
years and bolometric amplitudes of 1-2 mag (e.g., Blom- 
maert 1992). The bolometric luminosity variability affects 
the infrared properties through changes of the dust conden- 
sation radius and overall optical depth. Combining equation 
(0) with the luminosity relation L, — Anr^aT^ produces 

^ (lifer (^j-- 

A similar relation for galactic nuclei was noted in the recent 
numerical work of Laor & Draine (1993). Luminosity vari- 
ability induces movement of the dust condensation point: 
the shell's inner boundary is furthest from the star at maxi- 
mum light. This variation in linear scale-size of the envelope 
translates directly to its angular scale size, 9i. Envelopes are 
larger during maximum luminosity. It is noteworthy that 
variations of stellar radius and temperature during the stel- 
lar cycle do not enter independently, only through the overall 
luminosity. 

Such dependence of envelope size on phase of the light- 
curve has been detected. McCarthy, Howell & Low (1978) 
found that Mira (o Ceti) was partially resolved at 10.2 ^.m 
near maximum visual phase but unresolved at minimum 
phase. McCarthy, Howell & Low (1980) measured the size of 
IRC-l-10216 at 5 \im and obtained larger value at maximum 
than at minimum of the light-curve. Danchi et al. (1990, 
1994) reobserved this source at 11 ^.m and confirmed these 
results. In a separate study (Ivezic & Elitzur 1996; hereafter 
IE96) we analyse in detail all observations of IRC-l-10216 
and show that the angular diameter of its dust condensation 
zone varies between ~ 0.3 and ~ 0.5 arcsec as the luminosity 
changes from minimum to maximum. 

Variations of ri during the stellar cycle affect the dust 
optical depth so that r should be smallest at maximum light, 
largest at minimum. If the grain optical properties do not 
change during the stellar cycle, the overall optical depth is 
proportional to the dust column density. For constant mass- 
loss rate M and velocity i;i at the base of the envelope, the 
dust column density scales with r^^ . Together with equation 
(P), this implies r cx L^^^^. Since both M and vi may vary 
during the stellar cycle, the dependence of r on L, can be 
expected to be somewhat different but still show the same 
trend. From detailed analysis we find that, in IRC-l-10216, 
r oc Lz"'^ (IE96). Because of the variation of r, IR spec- 
tral properties should vary with the light curve, correspond- 
ing to redder radiation at minimum light (largest optical 
depth), bluer at maximum. Indeed, ever since the work of 
Pettit & Nicholson (1933), such behaviour has manifested 
itself in many observational correlations. In particular, Le 
Bertre (1988a,b) finds that in several dusty envelopes the 
wavelength of peak emission is shorter at maximum than 
at minimum light. In general, envelopes are bluer around 
maximum light. 

The phase dependence of optical depth can be anal- 
ysed quantitatively with the aid of photometric amplitudes 
A = 2.5 log(F""'"'/J'"^'"), where f and f""'" are fluxes at 
maximum and minimum light, respectively. At each wave- 
length, the photometric amplitude A\ is related to the am- 
plitude ylboi of the bolometric flux F via 



A\ = Ahoi + a\ , (7) 

where a\ is the amplitude of normalized flux fx — F\/F. 
Because of scaling, fx does not depend separately on lu- 
minosity, mass-loss rate, etc., only on overall optical depth 
(IE95). Therefore, ax is a unique function of r, specifically 

'^^W = 2-51og^^, (8) 

where r™^'' and r™'" are the optical depths at maximum 
and minimum light, respectively. That is, for a given dust 
chemical composition, ax too does not depend on individ- 
ual properties other than optical depth. Our detailed model 
calculations show that ax has a simple behaviour in two 
distinct limits. At wavelengths for which r < 1, typically 
A ^ 20 ^.m, ax oc — Aboi. At wavelengths for which emission 
can be neglected, typically A 3 \im (wavelengths shorter 
than the peak of the Planck function of the dust formation 
temperature Ti), oa oc ta. From these results we find that 

C Aboi + 0.7rA A ^ 3 ^m 

Ax ^ { (9) 
t (0.5 - 0.7) X Aboi A ^ 20 ^m. 

The photometric amplitudes vary linearly with r in the near- 
IR, becoming wavelength independent in the mid- and far- 
IR. Since Ahoi varies little among stars while r varies a lot, 
the range of photometric amplitudes can be expected to be 
large at short wavelengths, small at long ones. In addition, 
r decreases with wavelength in the near-IR, thus the ampli- 
tudes should decrease too. 

These simple results readily explain the available obser- 
vations. Short-wavelength amplitudes indeed decrease with 
wavelength and visual amplitudes are usually much larger 
than bolometric ones (e.g. Hetzler 1936; Lockwood & Wing 
1971). Furthermore, analysis of the Harvey et al. (1974) ob- 
servations shows that near-IR amplitudes of different stars 
increase with overall optical depth, and for a given star de- 
crease with wavelength. This behaviour is conflrmed by sub- 
sequent observations (e.g. DiCiacomo et al. 1991). Finally, 
the spread among stars is much larger for visual than for 
bolometric amplitudes (e.g. Forrest, Gillet & Stein 1975). 
For long wavelengths, Herman et al. (1984) flnd that the 
variation of OH maser luminosity, and presumably far-IR 
flux at 35 ^im, is about half that of the bolometric luminos- 
ity. 

Dust spectral features can display a more complex vari- 
ability with the stellar phase. Depending on initial optical 
depth, the strength of the 9.7-p.m silicate emission feature 
can either increase or decrease for the same change of r: in 
optically thin envelopes the strength increases with r un- 
til it reaches maximum at some optical depth Tm (~ 1), 
after which it decreases because of self-absorption (IE95). 
As a result, in envelopes with r < Tm the 9.7-|i.m emis- 
sion feature is strongest during minimum light (the optical 
depth is largest), with the reverse behaviour in envelopes 
with r > Tm. In envelopes with r ~ Tm, the emission fea- 
ture can remain almost unchanged during the stellar cycle. 
Little-Marenin, Staley & Stencel (1993) find examples for 
each type of behaviour in a number of sources that display 
the 9.7-jxm emission feature. On the other hand, when this 
feature is in absorption, its depth always increases with r 
(Scoville & Kwan 1976; IE95). Thus it is always deepest at 
minimum light, as observed (e.g. Engels et al. 1983). 
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4 OBSERVATIONAL IMPLICATIONS 

The source intensity is a unique function of various intrinsic 
properties such as luminosity, mass-loss rate, dust density, 
dust-to-gas ratio, etc. However, because of scaling, intrinsic 
properties do not enter separately, only indirectly through 
their effect in determining r and ri. For given grain optical 
properties, the surface brightness is a unique function of the 
scaled variable 6/ri, essentially determined by overall optical 
depth. Since optical depths involve only intrinsic properties, 
the same applies to the self-similar profiles displayed in Figs 
1 and 2. Contact with observations brings in the distance 
to the source, which fixes the scales of angular size 81 and 
bolometric flux F; that is, the distance sets the scale of con- 
version from linear to angular sizes and from intrinsic to 
observed fluxes. An additional dependence on the intrinsic 
parameters Ti and T, does exist, but it is rather weak and 
affects only optically thin envelopes at short (;$ 7 (i.m) wave- 
lengths. Furthermore, compared with the variation range of 
T and 01, Ti and T, are practically constant among late-type 
stars. As a result, the only non-radiative system properties 
that generally can be determined by the most detailed IR 
observations are r (including its wavelength dependence), a 
purely intrinsic property, and 9i , which sets the scales for ob- 
servations. We discuss now which observational techniques 
are most suitable for determining these two quantities. 

The least detailed observations are unresolved flux mea- 
surements, determining the spectral energy distribution Fx. 
This function is the product of the spectral shape f\ and 
the bolometric flux F. For given grains, the intrinsic prop- 
erty fx is controlled by r and is independent of 9i, thus it 
cannot be used to determine 61. However, for this very rea- 
son - its dependence on a single quantity - spectral shape 
analysis is the most reliable method for determining r and 
its wavelength dependence, resulting in the best handle on 
the dust chemical composition. The additional information 
contained in the bolometric flux allows an estimate of 9i, 
described below, albeit an indirect one. 

Direct determination of 9i requires spatially resolved 
observations. However, as shown in Section 2, 61 and r en- 
ter independently only in optically thin envelopes. In such 
sources, 61 can be directly read off the radial profiles from 
the position of the dust features, and r determined from 
the intensity ratio of stellar and dust components (cf. Fig. 
1). This determination of r is considerably less certain than 
that of 01. As a result, the surface brightness distribution is 
always a poor indicator of r, suggesting that the most prac- 
tical approach involves a two-step data analysis procedure. 
In the first step r is determined from a fit to the spectral 
shape, and in the second the resulting model prediction for 
the surface brightness is used to determine 61 . We have used 
such a two-step approach to analyse the IR observations of 
IRC-l-10216, and with a single model obtained successful fits 
to all spectral and spatially resolved observations (IE96). 

We proceed now to a more detailed discussion of the 
capabilities of different high-resolution techniques. 

4.1 Direct imaging 

Direct imaging is the best method for determining 9i, from 
the location of the features corresponding to dust conden- 
sation. Unfortunately, these features are visible only for a 



limited range of optical depths. At large r the features are 
smeared out as the profile evolves into a bell shape. At small 
values of r, the features are so weak compared with the stel- 
lar component that they become undetectable at the dy- 
namic range of most current telescopes; note, however, that 
on telescopes such as HST, Keck or Gemini, a dynamic range 
of 1000 is possible in images at spatial resolutions better 
than 0.1 arcsec. As a result, at every wavelength the features 
are a useful indicator only when the optical depth is of order 
unity. Thus the features generally cannot be used at either 
end of the spectrum: at short wavelengths 2-3 [im), the 
dust contribution is minute in comparison with the stellar 
component when r ~ 1 so the features are hard to resolve; 
and, at long wavelengths (Jj 20-30 \im), t never reaches 
unity for the majority of late-type stars (IE95). Therefore, 
the mid-infrared is the best spectral range for resolving the 
dust condensation features in surface brightness profiles. 

Table 1 lists the late-type stars we have identified as the 
best candidates for observations that can potentially resolve 
the dust formation features. They were selected as the stars 
with the largest expected values of 9i among IRAS sources 
whose 10 \im optical depth is in the range 0.1-2. In searching 
the IRAS catalogue we were careful to include only sources 
with the highest fiux qualities in all four wavebands and un- 
contaminated by cirrus emission (see IE95 for details of the 
selection process and removal of contaminated sources) . The 
optical depth was determined from the IRAS fluxes through 
extensions of the techniques of IE95, described in the Ap- 
pendix. Estimating the angular size 61 is more involved. Our 
starting point is the relation 

^--"■^^H io-Bwm-O [-rr) 

obtained by combining equation (^) with L» = AnD^F; note 
that when 9i and F are directly measured in observations, 
this equation can be inverted to determine the dust conden- 
sation temperature Ti. The Appendix describes a method 
for estimating the bolometric flux F from the IRAS fluxes, 
and for simplicity we assume Ti — 800 K for all stars. The 
resulting estimates for 9i axe expected to be accurate to 
within a factor of 2 or so. Table 1 presents the flve sources 
with the largest 9i that we identified in each of the follow- 
ing LRS classes: 10 (featureless spectra), 20 (silicate 9.7-^im 
feature in emission) and 40 (11.3-n.m in emission, mixture 
of amorphous carbon and SiC grains). 



4.2 Visibility observations 

All the radiative information available is contained in the 
surface brightness distribution. Therefore, direct imaging 
would be in principle the ultimate method to determine the 
system properties. However, finite dynamical range limita- 
tions preclude direct imaging of systems with very thin en- 
velopes, where the intensity of dust emission is negligible 
in comparison with the stellar component. The structure of 
such systems can still be probed by visibility observations 
for the following reason. As is evident from equation (^, 
visibility is a measure of fiux rather than intensity. And be- 
cause of the much larger area of the dust emission, its flux 
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Table 1. Best candidates for imaging; Late-type stars identified from their IRAS fluxes as the most promising 
candidates for potential resolution of the dust condensation zone in future observations. Quantities tabulated in 
columns (d)— (f) are model-based estimates. 



IRAS Name 


Other Name'''' 


LRSt'') 






no*"' 


01 (f) 


05027-2158 


T Lep 


15 


157 


1.36 


0.23 


0.03 


11294-6257 


na 


15 


5 


3.95 


1.76 


0.05 


12380-1-5607 


Y UMa 


15 


193 


0.89 


0.15 


0.04 


172974-1747 


V833 Her 


14 


559 


2.24 


0.38 


0.05 


20120-4433 


RZ Sgr 


16 


38 


1.91 


0.32 


0.02 


03507-1-1115 


IK Tau 


26 


4630 


0.12 


0.39 


0.22 


04566-1-5606 


TX Cam 


27 


1640 


0.05 


0.17 


0.22 


07209-2540 


VY CMa 


24 


9920 


0.27 


0.87 


0.31 


18050-2213 


VX Sgr 


26 


2740 


0.11 


0.37 


0.20 


23558-1-5106 


R Cas 


24 


1340 


0.06 


0.20 


0.17 


09116-2439 


na 


42 


737 


1.43 


0.22 


0.07 


15082-4808 


na 


42 


793 


1.41 


0.21 


0.08 


17049-2440 


na 


42 


793 


1.79 


0.27 


0.08 


18240-1-2326 


na 


42 


731 


1.73 


0.26 


0.07 


23320-1-4316 


LP And 


42 


959 


1.09 


0.16 


0.09 



(a) na: IRAS source without known association 

(b) IRAS LRS class 

(c) IRAS flux at 12 \mi (Jy) 

(d) Flux-averaged optical depth 

(e) Optical depth at 10 [om 

(f ) Angular size of the dust condensation zone (arcsec) 



can exceed that of the star even when the relation between 
intensities is the reverse 

At wavelengths shorter than ~ 2 \im and albedo ^0.1, 
the emission component of the dust radiation field is much 
weaker than either the scattered component or direct stellar 
contribution because it requires dust warmer than plausi- 
ble condensation temperatures. The radiation field is mostly 
composed of the attenuated stellar component and scattered 
light, and we find from our detailed solutions that under 
these circumstances 



K ^ exp(-rs, 



(11) 



where Tsca is the scattering optical depth. This relation has a 
simple physical explanation and it is valid for Tsca ^ 2. With 
this relation, visibility observations can determine the wave- 
length dependence of Tsca, a variation that can then be used 
to estimate the grain size a. The scattering optical depth 
is constant at short wavelengths, changing its behaviour to 
Tsca oc A~* for A ^ Aaca = 27ta. Determining Aaca from the 
change in behaviour of Tsca yields an estimate for a. We 
have applied this method to the visibility observations of 
IRC-l-10216 and obtained a ~ 0.2 fi.m, in agreement with 
an independent estimate based on spectral shape analysis 
(IE96). 



* This also implies that the stellar component can dominate the 
intensity profile even when it is negligible in the visibility obser- 
vations. This is contrary to a conclusion by Griffin (1990). 



5 DISCUSSION 

As shown in our previous study (IE95), the dynamics and 
spectral energy distributions of dusty winds around late- 
type stars are controlled by the dust chemical composition 
and overall optical depth. Here we have shown that a com- 
plete description of observations requires one additional pa- 
rameter - the angular size of the dust condensation zone. 
Since the spectral shape is entirely independent of 6i it is 
the most reliable indicator of r, whereas spatially resolved 
observations provide the most useful measurement of the 
dust condensation zone. Thus the most practical method of 
analysis is to determine the optical depth from model fits to 
the spectral shape, and 6\ from high-resolution observations 
(either direct imaging or visibility). Because of scaling, the 
model results presented here and in IE95 provide sufficient 
coverage of parameter space for the analysis of virtually any 
late-type star. An example of such complete analysis is pre- 
sented in IE96 for IRC-H0216. 

Our models do not include the effects of hght travel time 
between the central star and the heated dust. The phase lags 
that such travel introduces between variations of the dust 
and stellar components are insignificant at short- and mid- 
infrared wavelengths, whose emission originates close to the 
dust condensation zone. At long wavelengths, however, the 
dust emission originates sufficiently far from the star that 
the phase lags can become appreciable. Analysis of this effect 
has been recently reported by Wright & Baganoff (1995), 
who show that it can be used to determine the distance to 
the star. It is important to note that these effects do not 
invalidate our steady-state modelling, as discussed in IE95. 

In optically thin envelopes, our model images have a 
bright rim at the dust condensation radius, leading to the 
only features, except for the stellar spike, in the surface 
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brightness distribution. Images of optically thick envelopes 
do not display any details. We did not consider here the 
interaction of the wind with the surrounding interstellar 
medium, which may produce a dense shell far from the star 
(IE95). Such a shell would produce additional features in the 
surface brightness profile at long wavelengths (<; 50 i^m), 
which may have been observed for some IRAS stars at 60 
and 100 [im (Young, Phillips & Knapp 1993). Indeed, from 
the observed angular sizes of these features Young et al. con- 
clude that they probably correspond to the interface of the 
wind and the interstellar medium. 
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Table Al. Constants for equation (Al). 



A (|Jm) 


Sil.W 


Oli.(b) 


am.Cf'^) 


SiC(d) 


2.2 


0.4 


0.4 


1.20 


1.10 
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3.2 


3.6 


0.17 


0.15 


12.0 


1.5 


1.9 


0.12 


0.18 


25.0 


1.2 


1.5 


0.05 


0.04 


60.0 


0.2 


0.2 


0.01 


0.01 



(a) Astronomical silicate 

(b) Mixture of 80 per cent (by mass) astronomical silicate and 20 
per cent crystalline olivine 

(c) Amorphous carbon 

(d) Mixture of 80 per cent (by mass) amorphous carbon and 20 
per cent SiC 
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APPENDIX A: OPTICAL DEPTHS AND 
BOLOMETRIC FLUX FROM IRAS DATA 

Dust emission from late-type stars with envelopes is con- 
trolled by a single parameter, the flux averaged optical depth 
rp. The dependence of various spectral quantities on tf is 
listed by IE95, who also describe a method for determining 
rp from /7?^5 fluxes. Here we provide some additional useful 
relations, derived from our full model calculations. 

Once rp is known, the optical depth at any wavelength 
can be obtained from 

Tx = txr^. (Al) 

The power n is 1 for all dust grains when rp ^ 1. When 
rp > 1, n is 0.5 for silicate-based grains and 1.7 for carbona- 
ceous grains. The constants t\ at selected wavelengths are 
listed in Table Al for various grain compositions. For other 
wavelengths, tx can be interpolated using fig. 1 of IE95. Ex- 
trapolation toward shorter wavelengths is risky because of 
the increased contribution of scattering. 

An additional useful correlation involves the bolometric 
flux and IRAS fluxes: 

F = C X lO-'Veo f |i^V Wm-^ (A2) 

Here F\2, F25 and Feo are IRAS fluxes in Jy, and C and 7 
are constants that depend on the dust chemical composition. 
For silicate-based grains C = 7 and 7 = 3. For carbonaceous 
grains C = 0.08 and 7 = 6 when F12 > 2F25 and C = 1.8, 
7 = 1.5 otherwise. 
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